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“Appropriate  measures  and  consistent  standard  for  high  energy  laser  beam 

quality” 

Dr.  T.  Sean  Ross  (AFRL/DELO)^ 

Dr.  William  P.  Latham  (AFRL/DELS)'’ 

3550  Aberdeen  Ave  SE 
Kirtland  AFB,  NM  87117 

Abstract 

Along  with  power  output  of  the  laser  system,  laser  optieal  quality  or  beam  quality  provides  a 
suitable  measure  of  performanee.  Power  and  beam  quality  are  standards  for  the  eomparison  of 
laser  systems  with  eaeh  other  and  against  the  mission  requirements.  An  understanding  of  the 
meaning  of  beam  quality  is  neeessary  to  eompletely  define  the  lasers  performanee  eapability. 

The  eurrent  state  of  our  eommunity  ineludes  a  multitude  of  different  and  not  well  understood 
beam  quality  measures:  M  ,  Strehl  Ratio,  Brightness,  Power-in-the-bueket,  ‘times  diffraetion 
limited’  and  mode  eontent  determined  by  a  variety  of  beam  radius  measures:  half  widths,  seeond 
moment  radius,  widths  at  1/e  or  1/e  points,  width  of  primary  lobe  ete.  Another  eomplieation  is 
that  different  elements  of  the  eommunity  use  different  measures  to  evaluate  optieal  quality 
eharaeteristies. 

This  paper  will  examine  the  underlying  assumptions  behind  eommon  measures  of  beam  quality 
and  eompare  the  various  measures  as  they  relate  to  beams  from  lasers  employing  stable  resonant 
optieal  eavities.  We  show  how  the  mode  eomposition  of  a  beam  depends  upon  prior 
determination  of  beam  radius  and  how  the  term  ‘times  diffraetion  limited’  ean  mean  different 
things  depending  upon  the  method  used  to  measure  beam  radius.  We  show  the  ambiguities  that 
arise  between  eertain  elasses  of  beams  and  measures  of  beam  quality  and  advoeate  for  a  laser 
beam  quality  standard  that  relates  direetly  to  mission  requirements. 

KEYWORDS:  High-power  lasers.  Laser  beam  quality.  Laser  standards, 
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Symbols  and  abbreviations 


A, 
a 

B, 

P, 

e,C 

CCD 

CID 

D, 

d, 

5 

Av 

E 

e, 

f, 

FWl/e2M 

Act), 

H 

i 

l, 
i 

ISO 

/, 

LIDAR 

N 

k 

n, 

NBA 

NIR 

V 

m, 
m, 

M, 

MCA 

MIR 

MRC 

P, 

P, 

P, 

Pbi 

Pbi 


area 

spatial  limit  of  knife  edge  sean 
brightness 
beam  quality 
mode  eoeffieient 

charge-coupled  device.  Common  in  digital  cameras, 
charged  induced  device.  Common  in  digital  cameras, 
diameter 

measured  spot  size 

kronecker  delta  function. 

bandwidth 

electric  field 

natural  base,  2.7128... 

focal  length 

full  width  measured  at  l/e2  of  maximum 
azimuthal  cylindrical  angle 
rms  wavefront  distortion 
Hermite  polynomial 
summation  index 

irradiance  (radiometric  nomenclature)  and/or  intensity  (physics  nomenclature) 

normalized  irradiance 

International  Standards  Organization 

azimuthal  mode  index  for  Laguerre  Gaussian  functions 

laser  imaging,  detection  and  ranging 

number 

knife  edge  function 
wavelength 

X  index  for  Hermite  Gaussian  functions 

noise  equivalent  aperture 

near  infrared 

frequency 

measured 

y  index  for  Hermite  Gaussian  functions 
laser  mode  quality 
a  laser  beam  quality  measure 
mission  compliant  area 
mid-infrared 

mission  requirements  curve 
power 

power  fraction 

radial  mode  index  for  Laguerre  Gaussian  functions 

Ideal  power  in  the  bucket 

ideal  power  fraction  in  the  bucket 
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Pbm,  measured  power  in  the  bueket 

Pbm  measured  power  fraetion  in  the  bueket 

PIB  power  in  the  bueket 

n,  3.14159... 

r  radius 

rms  root  mean  square 

S,  Strehl  ratio 

s,  souree 

G  varianee 

t,  target 

0,  divergenee  angle  of  envelope 

0,  angle 

01  beam  divergenee  in  x  or  y  direetion 

0c  beam  divergenee  for  a  eireular  profile 

0g  beam  divergenee  of  a  eireular  gaussian  profile 

0g/  beam  divergenee  measured  at  the  lens  of  a  gaussian  profile 

0gi  beam  divergenee  in  the  x  or  y  direetion  of  a  gaussian  profile 

0g//  beam  divergenee  measured  at  the  lens  in  the  x  or  y  direetion  of  a  gaussian  profile 

0m  beam  divergenee  of  a  TEM00  gaussian  mode,  same  as  0g 

0h  beam  divergenee  of  a  higher  order  gaussian  mode 

TEM,  transverse  eleetro-magnetie 

TEMnm  the  nth,  mth  TEM  mode 

u  orthonormal  basis  funetion 

W,  measured  mode  radius 

w,  beam  radius 

Wo,  beam  waist 

Q  solid  angle 

XDE  times  the  diffraetion  limit 

x, y  transverse  distanees 

z,  axial  distanee 

zr,  Rayleigh  range 


Background 

The  output  of  a  laser  is  nearly  monoehromatie  and  extremely  eoherent.  Narrow  linewidth  or 
long  eoherenee  length  is  the  primary  operational  oharaeteristie  of  a  laser  and  it  allows  the  laser 
output  beam  to  be  foeused  to  a  tiny  spot.  Eor  a  thorough  diseussion  of  formal  eoherenee  theory 
and  relationships  see  refs  [1]  and  [2].  The  optieal  quality  or  beam  quality  of  a  laser  is  a  measure 
of  the  laser’s  foeusability.  More  importantly,  if  the  optieal  quality  of  the  laser  is  exeellent,  the 
laser  beam  is  eontrollable  and  understandable.  Many  solid-state  lasers  employ  a  resonant  eavity 
that  produees  a  beam  profile  made  up  of  a  series  of  Hermite-  or  Eagurre-Gaussian  modes.  Eor 
the  ease  of  a  stable  optieal  resonator,  the  laser  foeusability  is  eompared  to  the  foeusability  of  a 
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monoehromatie  field  with  the  lowest  loss  Gaussian  beam  spatial  profile  as  an  ideal  standard  to 
determine  its  optieal  quality.  If  the  eleetromagnetie  field  representing  the  modes  of  the  laser  is 
exaetly  known,  the  shape  and  power  in  the  beam  eould  be  evaluated  anywhere  within  the  foeal 
volume  of  that  beam.  However,  all  of  the  parameters  of  the  beam  are  not  generally  measured. 
Moreover,  a  solid-state  laser  generally  ineludes  thermally  indueed  self-foeusing  effeets  within 
the  laser  gain  material  as  the  laser  heats  up.  Although  there  are  analytie  ways  to  estimate  these 
effeets,  it  is  important  for  higher  laser  powers  to  make  an  empirieal  measure  of  the  solid-state 
laser’s  optieal  quality  under  operating  eonditions  to  evaluate  the  aetual  as-built  laser 
performanee.  The  two  eritieal  performanee  parameters  of  a  solid-state  laser  that  must  be 
determined  as  the  direeted  energy  eommunity  develops  solid-state  lasers  with  higher  powers  are 
laser  output  power  and  the  laser  optieal  quality.  The  purpose  of  this  paper  is  to  explain  the 
meaning  of  optieal  quality  and  mode  quality  in  basie  eoneeptual  terms,  so  that  it  is 
understandable  to  a  broad  speetrum  of  the  laser  development  eommunity  and  to  point  out  the 
assumptions,  strengths  and  weaknesses  of  various  measures  of  laser  beam  quality  to  enable  an 
informed  ehoiee  of  the  mission-appropriate  measure. 

Measures  of  Optical  Quality 

Optieal  quality  affeets  a  laser  beam’s  foeusability.  Measures  of  optieal  quality  are  divided  into 
two  eategories.  The  first  eategory  is  empirieal  or  measurable  quantities,  sueh  as  total  laser 
output  power  or  energy,  laser  linewidth,  foeal  spot  size,  far  field  peak  irradianee,  and  eneireled 
power  or  energy  in  the  foeal  spot,  whieh  is  sometimes  ealled  the  power  in  the  bueket.  Also 
ineluded  in  this  eategory  are  parameters,  whieh  are  direetly  ealeulated  using  these  measured 
values,  sueh  as  beam  divergenee,  eoherenee  length,  eoherenee  width  or  area,  and  brightness. 
These  parameters  quantify  the  lasers  performanee  empirieally  without  eomparison  to  a  standard. 
A  seeond  eategory  involves  relative  parameters,  whieh  eompare  the  foeusability  of  an  aetual 
laser  beam  to  the  foeusability  of  an  ideal  standard  laser  beam,  sueh  as  M  ,  beam  quality,  mode 
quality,  and  Strehl  ratio^. 

Beam  Divergence 

Diffraetion  is  the  name  given  to  the  angular  spreading  or  divergenee  of  light.  A  measure  of  the 
diffraetive  spread  of  a  laser  beam  is  the  ratio  of  the  average  wavelength  to  the  half  size  of  the 
laser  beam.  Various  expressions  for  the  beam  divergenee  angles^  are  listed  in  Table  1. 
Experimentally,  the  foeused  spot  diameter  or  spot  width  ean  be  measured  for  a  given  laser.  An 
empirieal  beam  divergenee  is  defined  to  be  the  ratio  of  the  spot  diameter  or  spot  width  to  the 
distanee  between  the  foeal  plane  and  lens. 


Table  1  Beam  divergence  angles 


Expression 

Deseription 

Beam  Divergenee: 

Divergenee  angle  in  the  Tth  direetion,  I  =  x  or  y,  for  the  fundamental  mode 

er2>./Di 

(TEMoo  mode)  due  to  a  plane  wavefront  ineident  upon  a  reetangular  aperture  of 
dimension  Dx  x  Dy.  About  81  %  of  the  total  energy  is  eontained  in  a  reetangular 
spot  defined  by  these  divergenee  angles. 

Q,=  2A4XfD 

Divergenee  angle  for  TEMqo  mode  due  to  a  plane  wavefront  ineident  upon  a 
eireular  aperture  of  diameter  D.  About  84.5  %  of  the  total  energy  is  eontained  in  a 
eireular  spot  defined  by  this  divergenee  angle. 

0g=  D/f  =  2V(71Wo) 

Divergenee  angle  for  TEMqo  mode  in  the  far- field  of  a  eireular  Gaussian  beam  of 
waist  radius  Wq.  About  86.5  %  of  the  total  energy  is  eontained  in  a  spot  defined 
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0g;=  2wo/f  =  4>,/(7iD,) 

by  this  divergenee  angle.  0g  is  measured  at  the  beam  waist. 

Divergenee  angle  for  TEMqo  mode  of  a  eireular  Gaussian  beam.  0g;  is  measured  at 

Bgi^  D;,/f  =  2>./(7iWoi) 

the  foeusing  lens. 

Divergenee  angle  in  the  Tth  direetion,  i  =  x  or  y,  for  TEMqo  mode  in  the  far-field 

0g/i=  2wo/f  =  2>^/(7iD/,) 

of  a  reetangular  Gaussian  beam.  0gi  is  measured  at  the  beam  waist. 

Divergenee  angle  in  the  Tth  direetion,  i  =  x  or  y,  for  TEMqo  mode  of  a  reetangular 

0m=d/f 

Gaussian  beam.  0g;,  is  measured  at  the  foeusing  lens. 

Divergenee  angle  for  any  laser  beam,  ds  is  the  measured  spot  size  and  f  is  the 

for  0g  (see  above) 

foeal  length  of  the  foeusing  lens. 

Divergenee  of  the  TEMoo  Gaussian  mode  that  is  the  fundamental 

0h  =  M  0, 

mode  for  the  higher  order  mode  of  interest. 

Divergenee  of  the  higher  order  Gaussian  mode. 

Strehl  Ratio 

Strehl  Ratio"^,  is  most  eommonly  used  in  the  astronomy  eommunity.  Astronomers  often  image 
point  objeets  loeated  in  the  extreme  far  field.  Strehl  is  a  ratio  between  the  peak  irradianee  of  a 
measured  signal  and  the  ealeulated  peak  irradianee  of  an  aberration  free  signal.  In  the  laser 
eommunity,  Strehl  eustomarily  means  the  ratio  between  the  peak  irradianee  of  a  laser  beam  and 
the  peak  irradianee  of  a  zero  order  Gaussian  or  other  ideal  beam  with  the  same  power  and  beam 
radius  as  in  equation  (1). 


ra?^x~ideal 

(1) 

The  Strehl  Ratio  is  eommonly  used  in  the  analysis  and  modelling  of  optieal  system  performanee. 
In  the  presenee  of  optieal  aberrations,  mirror  distortions,  gain  medium  imperfeetions,  and,  or 
beam  jitter,  the  far  field  peak  irradianee  is  degraded.  If  these  optieal  distortions  are  small,  the 
ratio  of  the  aetual  or  aberrated  far  field  peak  irradianee,  1^^,  to  the  ideal  or  unaberrated  peak 
irradianee,  Iq,  is  the  Strehl  Ratio,  S,  given  in  Table  2.  The  Strehl  Ratio  has  the  advantage  that  it 
is  simple  and  ean  be  useful  in  troubleshooting  laser  optieal  performanee  problems  or  during  laser 
design  when  a  partieular  distortion  will  potentially  be  introdueed  by  the  system.  An  empirieal 
form  of  the  Strehl  ratio  is  given  by  measuring  the  aetual  far  field  peak  irradianee  and  ealeulating 
the  ideal  peak  irradianee  from  the  laser  output  power. 

Brightness 

There  are  many  definitions^’^  of  brightness  within  the  laser  and  opties  eommunity.  In  general,  the 
brightness  (B)  for  a  eoherent  laser  souree  of  output  area,  As=Asource,  and  total  power,  P,  is  given 
by  B  =  P/  (As  G)  where  Q  is  the  solid  angle  subtended  by  the  area  in  the  far  field  or  observation 
plane  a  distanee  z  =  f  (f=l  foeal  length)  from  the  laser  souree  as  in  equation  (2).  The  solid  angle 

IS  given  by  O  ^target  AtargeJ^  ^  /As. 


5  = 


Z^P 


f2p 


/to  A  O 

-^source^  ^ toxget 


A  A  A  A 

-^source  arg  et 


(2) 
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There  is  no  ambiguity  in  ^source-  It  is  defined  in  terms  of  the  hard  aperture.  This  is  not  the  ease 
with  Atarget  whieh  ean  be  defined  in  terms  of  seeond  moment,  first  lobe  or  other  beam  radius 
measure.  With  suitable  eare.  Brightness  ean  be  a  eonsistent  measure.  If  ones  mission  is  to 
effieiently  put  power  in  a  solid  angle,  sueh  as  for  LIDAR,  traeking  or  oommunieations  missions, 
then  Brightness  is  the  mission  appropriate  measure.  It  is  also  eustomary  to  make  Brightness  a 
relative  measure  by  removing  power,  in  whieh  ease  it  is  ealled  etendue.  One  advantage  to 
Brightness  and  etendue  is  that  they  are  unaffeeted  by  the  opties  of  a  system.  In  the  absenee  of 
aberrations,  absorption  and  turbulenee.  Brightness  is  a  eonserved  property  of  an  optieal  system. 
Various  expressions  for  brightness  are  given  in  Table  2,  and  a  diseussion  of  this  topie  is 
presented  in  ref  [3]. 

2 

Beam  Quality,  M  ,  or  Mode  Quality 

The  primary  measurable  quantity  that  determines  laser  foeusability  performanee  is  the  eneireled 
energy  or  power  in  a  small  region  around  the  foeal  spot.  Beams  propagating  from  hard  apertures 
have  a  well-defined  eentral  lobe  whieh  determines  their  foeal  spot.  Soft  apertured  beams,  sueh  as 
Gaussian  modes,  do  not  have  a  well-defined  foeal  spot,  so  measures  sueh  as  seeond  moment  or 
1/e  must  be  used.  If  the  power  or  energy  measured  in  the  small  area  around  the  foeal  spot  is  Ps 
and  the  total  laser  output  power  or  energy  is  P,  the  fraetion  of  power  or  energy  delivered  by  the 
laser  to  the  foeal  spot  is  p  =  Ps  /  P.  The  fraetional  power,  p,  or  energy  in  the  far  field  foeal  spot 
is  a  readily  measurable  and  meaningful  quantity.  If  the  laser  is  operating  in  a  single  lowest  loss 
mode,  the  fraetional  power  ean  be  near  the  maximum  obtainable  for  some  ideal  standard 
waveform.  For  a  uniformly  illuminated  eireular  aperture,  p  »  0.84.  For  a  reetangular  beam 
geometry,  p  «  0.81.  For  the  lowest  order  Gaussian  beam,  p  «  0.86  by  the  1/e  eriteria.  The 
beam  quality,  P  is  defined  to  be  the  square  root  of  the  ratio  of  the  fraetional  power  in  the  far  field 
spot  for  an  ideal  standard  beam  to  the  fraetional  power  in  the  aetual  laser  beam,  that  is  P  = 

1/2 

(Pideal/Paet)  •  P  =  1  if  the  aetual  beam  has  the  same  far  field  power  as  the  ideal  standard.  The 
physieal  meehanism  that  determines  the  mode  strueture  of  the  laser  beam  and  the  foeused  beam 
spread  is  diffraetion.  For  p  =  1,  the  laser  is  said  to  be  diffraetion-limited.  The  Strehl  Ratio  beam 

1/2 

quality,  Ps,  is  given  by  Ps  =  (1/S)  .  The  bare  eavity  modes  of  a  stable  resonant  eavity  are 

Gaussian  beams.  For  a  single  higher  order  Gaussian  mode,  the  beam  waist  is  M  times  the  beam 
waist  of  the  lowest  loss  Gaussian  beam,  the  TEMoo  mode.  The  far  field  peak  irradianee  and  the 

2 

power  within  a  small  area  around  the  foeal  spot  is  redueed  by  a  faetor  of  M  ,  that  is  the 

2  2  1/2  2 

Gaussian  beam  quality,  Pg  ,is  related  to  M  by  Pg  =  (M  )  =  M.  M  or  Pg  provides  a 

meaningful  measure  of  the  optieal  quality  when  the  beam  strueture  eonsists  of  some  eombination 
of  the  Gaussian  modes^  ^  The  produet  of  beam  waist  times  far  field  divergenee  is  known  as  the 
beam  parameter  produet  and  is  another  measure  of  optieal  quality  whieh  ean  be  related  to  M 
and  is  most  eommonly  used  to  measure  the  beam  quality  of  semieonduetor  lasers.  If  there  are 
several  sourees  of  aberration  or  beam  degradation,  the  total  beam  quality  is  equal  to  a  produet  of 
all  individual  eontributions,  so  that  P  =  fl  Pi. 

The  various  measures  of  determining  beam  radius  will  be  diseussed  later,  but  it  is  important  to 
mention  here  that  most  measures  of  beam  quality  are  highly  dependent  upon  the  method  ehosen 
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to  determine  beam  radius  and  that  there  ean  be  eonsiderable  variation  in  measured  beam  quality 
depending  upon  beam  radius  definition. 


Table  2  A  short  list  of  measures  of  optical  quality 


Optical  Quality  Parameter 

Expression 

M  and  Factors: 

M  =  Bh/Bg 

M"  =  Bhf  Bg,  Bhy/Bgy 

Gaussian  Beam  Quality 

Beam  Parameter  Product 

M  =  (2p+/+l)'^^  for  circular  TEMp/  mode 

Mx  =  (2m+l)*^^  in  the  x  direction  for  rectangular  TEMmn  mode 

My  =  (2n+l)*^^  in  the  y  direction  for  rectangular  TEMmn  mode 

Pg  *  M 

Wo  Q 

Brightness 

Total  Brightness 

Spectral  Brightness 

B  =  P/(A,Q)  = 

Bt  =  P/O  =  PAA^ 

B^  =  P/(OAv)  and  B;,  =  P/(A,  Q  Ek) 

Strehl  Ratio 

S.  R.  Beam  Quality 

S  =  U/Io  ~  l-(27iA(|)/>.)2~Exp[-(27iA(|)/U^] 

Ps  =  1/  s'®  ~  1+(7iA(|)/>^)^  ~  Exp[+(7iA(|)/>^)^] 

Encircled  Power  Ratio  in  a 
Circle  of  Radius  r^ 

P  ~  (Pideal/pactual) 

Total  Beam  Quality 

3,  =  n,  B, 

Reduced  Brightness 

Bg  =  (l/3^)(PAUU)  =  f  h 

Power  in  the  bucket  (PIB) 

Power  in  the  bueket  is  simply  the  sum  of  the  total  power  within  a  partieular  area.  Usually  the 
area  is  eireular.  There  are  two  basie  means  for  determining  the  size  of  the  eirele.  The  first  is 
simply  to  use  the  aetual  size  appropriate  to  a  mission  target.  If  this  is  not  known  or  may  be 
variable,  then  one  may  use  a  size  based  on  the  diffraetion  properties  of  the  laser  output  aperture. 
The  radius  of  the  bueket  is  ealled  rb.  A  diseussion  of  how  rb  might  be  ehosen  for  a  partieular 
definition  of  beam  quality  is  given  in  a  later  seetion.  The  normalized  power  in  the  bueket  ean  be 
ealeulated  as  the  fraetion  of  output  power  that  ends  up  inside  the  target  eirele. 


r  7.71 

j  j  ^actual  (r,  (jfrdrdcj) 

NPIB  = 

j  \  Actual  {r,  (jfrdrdcj) 

0  0 

(3) 

Note  that  the  normalized  power  in  the  bueket  for  some  “ideal”  beam  is  given  by  equation  (3) 
with  aetual  replaeed  by  ideal.  If  the  ideal  beam  is  ehosen  to  have  the  same  total  power  as  the 
aetual  beam,  the  ratio  eomparing  the  power  inside  the  target  eirele  to  an  ideal  beam  is  given  by 
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PIBR 


’■hlTT 

J  actual  {r,  (j))rdrd(j) 

0  0 _ 

'i2;r 

J  \  fdeal  (r,  f)rdrd^ 
0  0 


P  actual 


Pideal 


1 


(4) 


The  ratios  are  useful  for  eomparing  the  effieieney  of  different  systems  but  may  lead  to 
disagreements  over  the  basis  of  eomparison.  Promoters  of  stable  resonators  will  likely  want  a 
PIB  based  on  an  ideal  low  order  Gaussian  beam  whieh  raises  all  the  ambiguities  in  determining 
the  charaeteristic  radius  of  an  appropriate  comparison  beam.  Those  working  with  unstable 
resonators  will  likely  want  a  comparison  with  a  flat-top.  The  actual  power  is,  of  course,  what 
ends  up  accomplishing  the  mission  and  represents  the  bottom  line.  In  the  comparisons  that 
follow,  we  use  equation  (4)  with  a  target  circle  equal  to  three  times  the  second  moment  radius. 

Comparison  of  beam  quality  measures 

Figure  1,  Figure  2  and  Figure  3  show  comparisons  of  M  ,  Strehl,  Brightness  and  PIBR  for 
various  laser  beams  to  illustrate  first,  that  these  measures  of  beam  quality  are  distinct  and  ought 
not  to  be  conflated  in  our  thinking  and  second,  to  provide  some  insight  to  the  researcher  to 
decide  which  measure  is  appropriate  for  a  given  application.  For  the  stable  resonator  modes  in 
Figure  1  and  Figure  2,  the  measures  of  beam  quality  were  calculated  using  the  z=0  plane  and  the 
z=3000m  plane  for  near  and  far  field.  No  atmospheric  distortions  were  included.  The  beam 
itself  was  chosen  to  have  a  z=0  (near  field)  second  moment  radius  of  3  cm.  M  was  calculated 
from  the  mode  coefficients  using  equation  (26).  M  is  always  greater  or  equal  to  1  and  will  be 
shown  with  numbers  on  the  right  hand  scale.  Strehl,  Brightness  and  PIBR  were  scaled  to  make 
their  maximum  value  1  and  are  shown  with  numbers  on  the  left  hand  scale.  Strehl  was 
calculated  from  the  peak  irradiances  of  near  and  far  field  beams.  The  Brightness  used  is  a 
relative  brightness,  sometimes  called  etendue,  based  on  the  physical  aperture  size  in  the  near 
field  and  the  second  moment  beam  radius  in  the  far  field.  This  was  so  the  brightness  would  have 
a  maximum  value  of  1  and  would  fit  on  the  same  graph  as  the  other  measures.  Relative  Power- 
in-the-Bucket  (PIBR)  was  calculated  using  the  total  power  inside  a  circle  with  radius  equal  to 
three  times  the  second  moment  waist  of  the  zero  order  Gaussian  beams  using  equation  (4).  This 
was  done  to  provide  a  fair  comparison  since  the  zero  order  Gaussian  beams  in  the  rest  of  the 
study  were  assumed  to  come  from  a  hard  aperture  with  radius  three  times  their  second  moment 
waist. 
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Figure  1,  Comparison  of  beam  quality  measures  for  a  beam  composed  of  TEMqo  and  TEMn  modes 


0  0.2  0.4  0.6  0.8  1 

Fraction  donut  mode 


Figure  2,  Comparison  of  beam  quality  measures  for  a  beam  composed  of 
TEMoo  and  TEMoi+TEMio  (donut)  mode 

Figure  1  and  Figure  2  show  M  ,  Strehl  ratio,  relative  brightness  and  PIB  as  a  function  of  percent 
of  higher  order  modes.  The  left  hand  side  of  these  figures  are  for  pure  zero  order  Gaussian 
modes,  TEMoo,  and  the  right  hand  sides  are  for  pure  order,  TEMn,  and  pure  TEMio  +TEM01, 
‘donut’  modes  as  shown  in  Eigure  11.  Eigure  1  shows  that  as  the  percentage  of  the  beam  that  is 
first  order  mode  increases  from  0  to  100%,  M  increases  from  1  to  3,  Strehl  drops  from  1  to  0, 
relative  Brightness  decreases  to  ~13%  while  relative  PIB  decreases  to  about  55%.  Eigure  2 

2 

shows  that  as  the  percentage  of  the  beam  that  is  ‘donut’  mode  increases  from  0  to  100%,  M 
increases  from  1  to  2,  Strehl  drops  from  1  to  0,  relative  Brightness  decreases  to  ~25%  while  PIB 
decreases  to  ~88%.  Comparing  two  systems  for  suitability  to  a  defense  mission,  is  a  donut  mode 
50%,  0%,  25%  or  88%  as  good  a  beam  as  a  zero  order  Gaussian?  Is  a  TEMn  beam  33%,  0%, 
13%  or  55%  as  good  a  beam  as  a  TEMoo  Gaussian?  The  differing  measures  of  beam  quality 
clearly  do  not  measure  the  same  thing. 
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Figure  3,  Comparison  of  beam  quality  measures  for  a  square  annulus  flat  top  beam  for  various  apertures 

Figure  3  shows  a  comparison  of  various  beam  quality  measures  vs.  the  outcoupling  size.  It  is 
important  to  emphasize  that  is  based  on  the  second  moment  measure  of  beam  radius. 
Theoretically,  second  moment  requires  an  infinite  plane  -  all  information  is  significant  and 
information  away  from  the  centroid  is  quadratically  more  and  more  significant.  here  is 
calculated  from  Fourier  propagation  on  a  1024x1024  numerical  array  with  spacing  of  0.2mm  -  a 
total  field  approximately  20cm  across.  The  far  field  second  moment  radius  calculated  is  a 
function  of  this  numerical  field,  similar  to  the  way  a  measured  M  measurement  will  be  a 
function  of  the  detector  size  and  noise  equivalent  aperture  as  will  be  discussed  in  the  section  on 
experimental  foundations  of  M  .  Differing  systems  will  measure  varying  M  values  for  the  same 
beam.  We  thus,  do  not  intend  to  convey  that  the  of  a  square  annulus  flat  top  is  28. 
Theoretically,  is  infinite  for  all  beams  from  hard  apertures.  We  do  intend  to  convey  that 
can  be  measured  to  be  28  in  conditions  which  match  those  of  our  Fourier  propagation  code.  In 
Figure  3,  increases  from  ~5  to  ~28  as  the  fractional  aperture  increases  from  0.05  to  0.95. 
Strehl  (compared  with  a  square  fiat  top  beam)  decreased  from  1  to  0.0,  relative  brightness 
decreased  from  0.25  to  zero  and  PIB  decreased  from  88%  to  15%.  In  comparing  a  square  fiat  top 
with  a  square  annulus  flat  top  with  a  fractional  aperture  of  0.4,  is  the  apertured  beam  70%,  75%, 
60%  or  97%  as  good  a  beam  as  the  square  flat  top? 

The  Foundations  of 

The  term  M  has  developed  out  of  the  theory  of  Gaussian  modes.  Gaussian  modes  are  the  bare 
cavity  solutions  for  stable  resonant  cavities.  Within  this  section,  the  Laguerre-  Gaussian  modes 
which  are  appropriate  for  radial  and  azimuthal  coordinates  will  be  discussed.  The  Cartesian 
coordinate  modes,  the  Hermite-Gaussian  modes,  are  given  below.  The  higher  order  radial 
Gaussian  modes  grow  in  size  as  a  function  of  the  radial  mode  index.  The  radial  size,  rn,  of  the 
nth  mode^°  is  approximated  by: 
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nxwQ,r^ 


=  Mw  = 


where  wo  is  the  radius  of  the  lowest  order,  O*’’  order,  Gaussian,  the  TEMqo  mode.  The  Laguerre- 
Gaussian  modes  and  Hermite-Gaussian  modes  are  given  in  referenee  [11].  We  have  ehosen  to 
label  the  sealing  faetor  for  the  mode  M  for  reasons  that  will  be  elear  in  the  next  seetion.  Some 
times  the  radial  mode  size  is  labeled  Wn.  The  lowest  order  Gaussian  mode  beam  radius  has  a 
quadratie  profile  as  a  funetion  of  axial  distanee  along  the  eenter  of  the  beam. 


W  =  w(z)  =  Wg 


1  + 


do{z-z,) 


.2  ^ 


Wn 


=  w, 


1  + 


,2  ^ 


J 


TW 


2 

0 


(6) 


where  wo  is  the  O*’’  order  Gaussian  beam  waist  radius,  6*0  is  the  0*  order  Gaussian  beam 
divergenee  angle  and  zr  is  the  Rayleigh  range.  The  higher  order  modes  follow  a  similar  quadratie 
profile  with  the  beam  radius  sealed  everywhere  in  z  by  the  sealing  faetor  M.  The  M  value  or 
the  mode  quality  is  determined  by  measuring  the  quadratie  envelope  of  a  partieular  aetual  laser 
beam  and  finding  the  fit  parameters  for  a  quadratie  eurve.  The  quadratie  envelope  ean  readily 
be  ealeulated  for  the  envelope  W{7).  The  expression  for  the  quadratie  that  follows  the  Gaussian 
beam  envelope  is 


W{z)  =  W, 


'0,{z-z^)Y 

I  K  J 


W  (z)  =  Mw[z^ 


(7) 

where  the  2"‘*  moment  Gaussian  beam  radius  for  the  0**'  order  Gaussian  is  w(z),  0  is  the 
divergenee  of  the  envelope,  and  the  beam  radius  for  the  aetual  beam,  W{z),  is  taken  to  be  the 
sealing  parameter,  M,  times  the  0*  order  Gaussian  as  given  in  equation  (5).  The  result  of  the  for 
the  envelope  quadratie  funetion  is 

A  tiWq 

(8) 

Note  that  we  have  ehosen  to  label  the  radius  for  the  0*  order  beam  with  lower  ease  symbols;  1) 
beam  radius  is  w(z)  and  2)  beam  divergenee  is  6*.  The  quadratie  envelope  parameters  for  the 
aetual  measured  beam  are  labeled  with  upper  ease  symbols;  1)  beam  radius  is  W{z)  and  2)  beam 
divergenee  is  0.  Determining  the  M  value  depends  on  what  partieular  measurement  seheme  is 
used  and  what  ehoiees  are  made  for  values  within  the  fitting  proeess.  Usually,  the  aetual  beam  is 
a  mixture  of  a  few  or  many  Gaussian  modes.  There  are  analytie  expressions  for  the  sealing 
faetor  for  single  higher  order  modes.  However,  there  is  not  a  unique  Gaussian  mode  series 
solution  for  a  given  eleetrie  field  profile  or  envelope.  In  the  absense  of  any  laser  fiuetuations  or 
noise,  three  measurements  of  the  beam  radius,  Wi,  and  beam  positions,  Zi,  at  three  loeations  are 
suffieient  to  determine  the  quadratie  envelope  of  the  beam  if  there  is  no  noise  or  beam 
fiuetuation  of  any  kind.  In  praetiee'"^,  many  measurements  must  be  taken.  If  the  beam  eonsists 
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of  a  single  higher  order  Gaussian  mode,  the  value  for  M  can  be  calculated  by  measuring  the  2"*^ 
moment  beam  radius  and  multiplying  by  the  mode  coefficients  as  in  equation  (26). 


Figure  4,  Determining  the  Parameters  for  the  Gaussian  Envelope 

2 

For  a  multimode  Gaussian  beam,  the  beam  also  has  a  quadratic  shape.  M  is  based  upon  the 
comparison  of  the  multimode  beam  to  an  ideal  TEMoo  beam.  However,  as  the  number  of  modes 
increases,  the  proper  beam  radius  of  that  ideal  TEMoo  beam  becomes  increasingly  hard  to 
determine  and  creates  inherent  uncertainty  in  the  measure  of  as  will  be  discussed  below.  The 
irradiance  for  the  one  dimensional  Hermite-Gaussian  modes,  TEMn  modes,  are  shown  in  Figure 
5. 
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Hermite  Gaussian  Modes 

Irradiance  -  Scaled  to  Laser  Power,  P 


wO  =  1  Mode  Index 
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Figure  5  One  Dimensional  Hermite-Gaussian  Mode  Irradiances 

The  formula  for  these  modes  is  given  in  equation  (21).  Note  that  the  modes  increase  in 
transverse  size  as  the  mode  index  increases  as  in  equation  (5).  After  normalizing  to  the  laser 
power,  P,  the  irradiance  of  the  lowest  order  Gaussian  mode  is  given  by 


I(r) 


2P 

TTW^ 


2p  2{x^+y^) 

_ p 

2  ^ 

TTW 


(9) 


where  the  beam  radius  as  a  function  of  the  longitudinal  coordinate,  z,  and  the  beam  divergence 
angle  are  given  by  equation  (6).  Since  the  higher  order  Gaussian  modes  increase  in  transverse 
size  with  increasing  index,  the  irradiance  envelope  for  a  higher  order  mode  or  multimode  beam 
can  be  approximated  by  the  lowest  loss  Gaussian  that  is  larger  in  radius.  By  normalizing  to  the 
same  laser  power,  the  larger  size  mode  envelope  is  given  by 


2P 

71  mW 


2d 

g  M^mP' 


2P 

- W  =  Mw 

TlW^ 


(10) 

where  the  beam  radius  and  beam  divergence  angle  are  given  in  equations  (7)  and  (8).  Once 
again,  the  expression  can  readily  be  converted  to  x  and  y  coordinates  by  using  r^=x^+y^.  For 
non-symmetric  cases,  the  envelope  irradiance  can  be  written  as  a  product  of  an  irradiance  in  x 
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and  an  irradianee  in  y  that  are  scaled  by  Mx  and  My.  In  that  case,  the  mode  quality  value,  ,  is  a 
product  as  follows 


M"  =  M^M^, 


(11) 

The  mode  quality  or  M  value  is  directly  related  to  the  quadratic  envelope  of  the  beam  as 
discussed  in  the  previous  section.  Another  measure  of  optical  quality  is  the  power  in  the  bucket 
as  given  in  equation  (3).  Here,  we  take  the  bucket  radius  to  be  the  1/e^  point  in  the  lowest  loss 
Gaussian  far  field  irradianee  pattern.  At  this  bucket  radius,  the  power  in  the  bucket  in  the  lowest 
loss  Gaussian  mode  is  given  by'^ 


P„  =  />  X  (1  -  A);  =  1  -  «  86.4% 


(12) 

That  is,  eighty-six  percent  of  the  total  laser  power  is  delivered  into  the  bucket  for  the  lowest  loss 
mode.  For  the  larger  radius  mode  envelope,  the  power  in  the  bucket  is  given  by 


f 


A  =  ^x 


\-e 


2  A 


lPbm=^-e 


(13) 

2  2 

Note  that  equation  (13)  reduces  to  equation  (12)  for  M  =1.  When  M  >1,  the  normalized  power 
in  the  bucket  is  less  than  eighty-six  percent.  Finally,  the  beam  quality  is  usually  taken  to  be  the 
square  root  of  the  ratio  of  the  ideal  standard  beam  power  in  the  bucket  to  the  actual  beam  power 
in  the  bucket  for  the  same  bucket  radius.  Thus,  the  beam  quality  of  the  envelope  beam  can  be 
given  by 


(14) 

The  transverse  envelope  function  has  been  used  to  determine  a  “beam  quality”  based  on  a 
comparison  to  the  power  in  the  bucket  of  the  lowest  loss  Gaussian  beam  as  an  ideal  standard. 

This  gives  a  definition  of  beam  quality  for  a  multimode  Gaussian  beam.  Asymmetric  beams  can 
readily  be  analyzed  in  a  similar  way. 

There  is  an  important  conceptual  subtlety  involved  in  the  use  of  equation  (13).  As  one  departs 
from  single  mode  content,  it  becomes  less  and  less  clear  what  the  proper  radius  of  the  lowest 
order  embedded  Gaussian  is.  As  one  progresses  to  fiat  top  shapes,  there  are  many  different, 
reasonable  criteria  on  which  to  make  this  comparison  and  none  of  them  are  totally  convincing 
because  a  fiat  top  has  a  hard  aperture  while  Gaussian  modes  are  all  of  infinite  extent.  By  suitable 
comparison,  a  fiat  top  can  look  terrible  or  very  good  in  comparison  to  a  TEMoo  Gaussian.  If  we 
judge  each  beam  against  itself,  rather  than  by  a  tenuous  comparison  to  an  ideal,  then  we  find  that 
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a  TEMoo  Gaussian  is  not  vastly  superior  to  other  beam  shapes.  For  example,  A  TEMoo  Gaussian 
has  86%  of  its  energy  inside  a  1/e^  bucket  while  flat  top  beams  have  approximately  77%  of  their 
energy  inside  this  radius,  yet  flat-tops  have  an  infinite  M^.  The  general  rule  for  the  use  of 
and  equation  (13)  is  that  they  only  apply  to  beams  which  can  be  fully  and  uniquely  characterized 
by  a  few  Gaussian  modes.  Reference  13,  for  example,  reports  that  is  a  unique  and  meaningful 
measure  for  stable  resonator  beams  with  non-annular  round  beams  up  to  M^<3.2  .  In  the 
following  sections,  we  discuss  some  of  the  experimental  and  conceptual  problems  associated 
with  the  use  of  M  . 

M^-What  your  beam  analyzer’s  manual  didn’t  tell  you 

M  is  probably  the  most  popular  measure  of  laser  beam  quality.  It  is  also  one  of  the  most 
experimentally  misused  and  inconsistently  applied  measures.  Many  scientists  have  their  own  pet 
methods  for  measuring  M  :  some  use  a  knife  edge;  some  take  one  measurement  at  focus  and  one 
in  the  far  field;  some  measure  the  transmission  through  a  single  mode  optical  fiber;  some  relate  it 
to  Strehl  ratio  or  other  beam  quality  measures.  Many  researchers  rely  on  commercial  “black¬ 
box”  devices  and  accept  the  manufacturers’  assurances  that  the  device  actually  measures  M^. 

One  cannot  rigorously  examine  the  methods  these  various  commercial  devices  use  because  they 
are  hidden  behind  the  word  “proprietary”.  One  must  trust  the  manufacturer’s  sales  literature.  A 
proper  M  measurement  is  difficult  to  take.  Fortunately,  there  is  an  international  standard 
which  specifies  how  the  measurement  is  to  be  taken.  A  proper  M^  measurement  takes:  the  2"‘^ 
moment  beam  radius  measured  in  at  least  10  places  in  the  far  field  and  through  the  focus,  then  fit 
to  an  ideal  Gaussian  equation  in  terms  of  Rayleigh  Range,  focus  location  and  M^.  Further,  the 
experimenter  must  be  able  to  understand  the  sources  of  uncertainty  in  the  measurement  to  be 
able  to  correctly  assign  error  values,  something  sadly  lacking  from  most  commercial  “black  box” 
beam  analyzers. 

T rade  Offs 

The  heart  of  a  beam  quality  measurement  system  is  the  camera.  Silicon  based  charge  coupled 
devide  (CCD)  or  charge  induced  device  (CID)  cameras  are  the  most  common,  typically  have 
pixel  sizes  <  10|um,  detect  adequately  in  the  visible  and  near  infrared  and  are  relatively 
inexpensive.  For  the  mid-infrared,  there  are  screened  cameras  (with  pixel  bleed)  or  pyroelectric 
CCD  cameras  with  large  pixels  (100|um  x  100|um).  These  cameras  can  necessitate  the  use  of 
very  long  travel  stages  due  to  low  resolution.  A  long-travel  stage  allows  a  loose  focus,  a  smaller 
variation  in  irradiance  as  the  detector  is  moved  through  the  beam,  less  variation  in  signal  to  noise 
ratio  and  less  stringent  requirements  on  the  pixel  size  of  the  camera.  These  benefits  come  at  a 
price.  Short  travel  precision  stages  (<10cm)  are  relatively  inexpensive.  As  the  travel  distance 
increases  beyond  10  cm,  the  price  increases  significantly.  These  short  travel  stages  give  focused 
spot  sizes  appropriate  for  silicon  based  CCD  cameras  for  the  visible  and  NIR  but  are  often 
inadequate  in  the  MIR  where  screened  cameras  are  subject  to  pixel  bleed  and  pyroelectric  CCD 
cameras  that  have  very  large  pixels. 

A  further  implication  of  stage  length  is  the  size  of  the  image  in  the  CCD  camera  and  the  number 
of  pixels  across  the  beam.  If,  for  example,  we  fill  the  CCD  aperture  in  the  far  field,  we  may  have 
300^  pixels  across  the  beam.  With  a  short  travel  stage,  we  may  have  to  focus  down  to  only  10^ 
pixels  across  the  beam  with  a  corresponding  decrease  in  signal  to  noise  ratio.  The  reliability  of 
an  M^  measurement  is  in  the  near  focus  and  far  field  measurements.  Poor  signal  to  noise  near 
focus  is  a  source  of  error  in  the  final  measured  value  of  the  beam  radius.  Figure  6  shows  a 
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sample  beam  measured  in  the  far  field  and  near  focus.  The  accuracy  of  the  beam  radius 
measurement  is  far  lower  near  focus  than  in  the  far  field  simply  due  to  several  orders  of 
magnitude  fewer  pixels  across  the  beam.  Note  the  change  of  transverse  scale  in  Figure  6 
showing  the  significant  portion  of  each  beam. 


Figure  6,  Near  and  far  field  beam  measurements 


Pseudo-average 

The  measurement  of  M  is  plagued  by  a  number  of  theoretical-experimental  disconnects. 
Theoretically,  is  an  instantaneous  concept.  A  given  wavefront  at  a  given  instant  can  be 
assigned  an  M  value.  Experimentally,  M  is  a  time-averaged  measurement.  10  measurements 
with  intervening  stage  movements  filter  changes,  and  software  aperture  changes  takes  time.  The 
ISO  standard  specifies  that  a  laser  must  be  warmed  up  for  at  least  1  hour  prior  to  measurement. 
The  hope  is  that  the  laser  will  become  stable  enough  so  that  measurements  taken  on  different 

2 

parts  of  the  wave  train  will  yield  an  average  value  approximately  equal  to  the  instantaneous  M 
over  the  entire  wave  train.  If  the  laser  fluctuates  in  power  or  mode  content  during  the  course  of 
the  measurement,  the  results  are  less  valid.  If  the  laser  remains  constant  during  measurement  but 
fluctuates  afterwards,  the  results  were  valid,  but  not  useful.  is  thus  inappropriate  for  single¬ 
shot  lasers,  those  with  short  run  times  or  with  even  slow  fluctuations  in  power  or  mode  content. 

It  is  also  inappropriate  if  the  measurement  is  taken  under  different  operating  conditions  than  the 
system  will  be  used  in. 

Second  Moment  Beam  Radius 

There  are  many  ways  to  define  or  measure  the  width  of  a  laser  beam.  M  is  defined  in  terms  of 
the  second-moment  of  the  irradiance.  The  second  moment  is  calculated  by  weighting  the 
measured  irradiance  by  the  square  of  the  distance  from  the  centroid  of  the  beam. 

X  -  xf  I{x,y)  dxdy 

^^I{x,y)dxdy 


(15) 

The  use  of  the  second  moment  to  measure  beam  radius  creates  another  theoretical-experimental 
disconnect.  Theoretically,  the  second  moment  is  very  sound  and  can  be  applied  to  a  wide  variety 
of  beam  shapes.  Experimentally,  it  is  problematic  because  small  amounts  of  noise  away  from 
the  beam  are  weighted  more  heavily  than  the  actual  signal  you  are  measuring.  There  are  a 
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number  of  steps  that  must  be  taken  as  a  eonsequenee  of  the  ehoiee  of  seeond  moment  definition 
of  beam  radius.  The  first  is  due  to  the  faet  that  CCD  eameras  only  return  positive  signals.  An 
artifieial  zero  must  be  ehosen  sueh  that  the  seeond  moment  of  the  noise  is  zero.  The  seeond  is 
that  no  extraneous  information  should  be  taken  into  ones  measurement.  Extra  information  only 
eontributes  to  error.  Determination  of  what  eonstitutes  “extra  information”  will  be  diseussed  in 
the  seetion  on  noise  equivalent  aperture. 

Another  implieation  of  the  seeond-moment  heavily  weighting  data  from  the  wings  of  a  beam  is 
that  annular  beams  sueh  as  those  from  unstable  resonators  are  measured  to  have  inordinately 
large  beam  radii  and  suffer  in  eomparison  with  other  lasers.  Many  high  power  laser  systems  use 
unstable  resonators.  Defining  beam  quality  speeifieations  in  terms  of  for  a  high  power 
unstable  resonator  may  give  the  laser  designer  an  impossible  task  that  will  not  signifieantly 
enhanee  the  ability  of  the  system  to  aeeomplish  its  intended  mission.  For  example,  take  the  ease 
of  a  beam  with  of  1 .0,  entirely  O*’’  order.  Compare  this  with  the  beam  of  M^=2  eomposed  of 
95%  TEMoo  mode  and  5%  TEM20+TEM20  mode  and  a  beam  of  M^=2  eomposed  of  75%  TEMoo 
mode  and  25%  TEMn,  as  shown  in  Figure  7.  The  beam  in  the  eenter  has  some  higher  order 
noise  that  will  result  in  5%  of  the  energy  diffraeting  away  from  the  intended  target.  The  beam  on 
the  right  also  has  an  of  2,  the  same  as  the  beam  in  the  eenter,  but  it  has  severe  impaet  on  the 
shape  of  the  foeal  spot  and  may  result  in  a  loss  of  up  to  25%  of  the  energy  of  the  laser.  This 
diserepaney  is  a  result  of  the  faet  that  is  fundamentally  linked  with  the  seeond  moment 
measure  of  beam  radius. 


Figure  7,  M^=l  beam  (left)  M^=2  beams  (center,  right)  in  tbe  far  field 


Noise  Equivalent  Aperture 

With  a  beam  analysis  paekage  of  suffieient  flexibility,  one  ean  experimentally  verily  that  the 
measured  seeond  moment  radius  is  a  strong  funetion  of  the  software  aperture  put  around  the  data 
returned  by  the  CCD  eamera.  This  is  due  to  the  faet  that  no  matter  how  earefully  one  sets  the 
artifieial  zero,  the  seeond  moment  of  the  noise  will  fiuetuate  slightly  around  zero.  The  non-zero 
seeond  moment  of  the  noise  will  signifieantly  alter  the  results  of  the  measurement.  In  the  results 
shown  in  Figure  8,  the  measured  seeond  moment  beam  radius  for  the  same  beam  from  a 
standard,  eommereial  Nd:YAG  laser  varies  from  0.45  mm  to  0.15  mm  depending  on  what 
software  aperture  is  ehosen,  plotted  as  a  fraetion  of  the  full-width  at  1/e^  maximum  beam  radius. 
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Beam  Radius  vs.  Numerical  Window 


Fraction  of  FWl/e2M  Window 

Figure  8,  Measurement  of  2"“*  moment  beam  radius  as  a  function  of  software  aperture 


The  discrete  nature  of  the  signals  returned  from  digital  cameras  has  some  long-reaching  effects 
on  the  measurement  of  beam  quality.  If  we  take  a  hypothetical  case,  shown  in  Figure  9  of  the 
column  averages  for  an  8  bit  camera,  the  signal  returned  will  fall  between  0  and  255.  Further 
assuming  that  the  appropriate  artificial  “zero”  is  at  a  pixel  value  of  50,  this  gives  a  maximum 
contrast  of  205;  1 .  A  pixel  value  of  5 1  is  barely  significant.  A  pixel  value  of  50  is  effectively 
“zero”.  Any  information  taken  beyond  the  point  at  which  the  signal  is  less  than  1/205*  of  the 
peak  is  extraneous  and  can  only  contain  noise.  This  defines  a  “noise  equivalent  aperture”, 
derived  in  equation  (16)  assuming  a  TEMoo  gaussian,  is  equal  to  ~2.3  beam  radii  in  our  example, 
which  must  be  electronically  placed  around  the  CCD  image  to  exclude  parts  which  can  only 
contain  noise.  This  raises  its  own  problem  in  that  one  must  know  the  beam  radius  and  the  mode 
content  in  order  to  calculate  the  noise  equivalent  aperture  in  order  to  measure  the  beam  radius 
prior  to  determining  M  which  gives  only  clues  as  to  the  mode  content.  In  practice,  this  means 
that  the  beam  radius  must  be  first  measured  by  a  non-second  moment  method  to  correctly  set  the 
noise  equivalent  aperture  prior  to  measuring  the  second  moment  beam  radius.  Unfortunately,  the 
fiat  portion  of  the  curve  in  Figure  8  and  the  correct  noise  equivalent  aperture  do  not  coincide.  If 
they  did,  then  small  errors  in  pre-determining  the  beam  radius  wouldn’t  matter.  The  correct 
noise  equivalent  aperture  is  on  a  sloped  portion  of  the  curve  and  thus  errors  in  the  pre¬ 
determination  of  the  beam  radius  result  in  errors  in  the  measured  second  moment  beam  radius. 
Note  that  NEA,  thus  defined,  is  the  largest  NEA  appropriate.  A  careful  characterization  of  the 
dark-current  noise  of  the  CCD  camera  will  show  an  rms  fluctuation  which  will  likely  be  greater 
than  the  1  pixel  value  used  here.  This  will  reduce  the  contrast  accordingly.  Contrast  =  (peak 
signal-artificial  zero)/  Max[rms  noise  fluctuation  of  camera,!]  .  Note  that  frame  averaging  can 
be  used,  but  does  not  have  quite  the  effect  that  might  be  desired  and  is  discussed  in  Appendix  A; 
Derivation  of  error  terms  in  measurement. 
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Smallestsignal 
1/  205  of  peak 


Figure  9  Image  contrast  for  a  hypothetical  CCD  camera 
Figure  9, 


2NEA^  1 

g  w"  = -  NEA  =  >^^111(205) 

contrast 

X  m - \  =  2.3  beam  radii 

NEA  =  w-^jmycontrast ) 

(16) 

Equation  (16)  is  strictly  correct  only  for  a  TEMoo  beam.  As  higher  order  modes  enter,  the  NEA 
must  be  re-evaluated  based  on  the  new  mode  content.  This  places  the  experimenter  once  again 
in  the  quandary  of  needing  to  know  what  he  is  attempting  to  measure  in  order  to  measure  it 
properly.  The  measurement  of  M  is  thus  limited  to  small  perturbations  on  the  basic  TEMoo 
profile.  The  higher  is  measured  to  be,  the  less  meaningful  it  becomes. 

Irradiance  dependent  beam  radius 

One  further  implication  of  the  available  contrast  in  a  CCD  image  is  that  the  noise  equivalent 
aperture  is  a  function  of  the  peak  signal.  If  the  peak  of  the  image  returns  a  pixel  value  of  200, 
the  noise  equivalent  aperture  will  be  different  than  if  the  peak  of  the  image  was  at  250  due  to  the 
change  in  contrast  of  the  image.  Eurther,  a  low  contrast  image  will  have  a  lower  signal  to  noise 
ratio  and  measure  the  beam  radius  with  greater  error  than  a  high  contrast  image.  This  has 
relevance  to  ones  filter  set.  Neutral  density  filters  are  commonly  used  to  attenuate  the  beam  to 
protect  the  CCD  camera.  The  ideal  case  would  be  a  continuous  range  of  attenuation  to  keep  the 
peak  signal  near  the  saturation  point  of  the  CCD  camera  for  all  measurements  to  maximize  the 
image  contrast  and  keep  the  noise  equivalent  aperture  constant.  Typically,  one  has  a  discrete  set 
of  filter  values.  The  discontinuities  in  neutral  density  attenuation  provide  another  source  of 
uncertainty  in  the  final  measured  value. 

Curve  fitting 

Once  the  series  of  beam  radius  measurements  has  been  taken,  the  data  must  be  fit  to  a  quadratic 
equation  in  terms  of  M^,  the  beam  waist  (wo)  and  the  focal  position  (zq).  Note  that  there  are  two 
forms  of  this  equation,  shown  as  equations  (17)  and  (18).  The  difference  is  in  how  wo,  the  beam 
waist,  is  interpreted.  If  Wo  is  interpreted  as  the  beam  radius  of  a  pure  Gaussian  TEMoo  mode 
embedded  in  the  beam,  then  equation  (17)  is  used.  If  wo  is  interpreted  as  the  smallest  measured 
beam  radius,  then  equation  (18)  is  used.  Since  the  measurement  of  M^  deals  with  measured,  not 
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theoretieal,  beam  radii,  equation  (18)  is  the  most  eommon  in  curve  fitting^^,  while  equation  (17) 
is  the  most  common  in  scientific  literature^^. 

w\z)  =  M‘'wl+[-^){z-zf 

(17) 

w\z)=wl+{f^){z-zj 

(18) 

The  ISO  standard  does  not  specify  a  numerical  method  but  does  recommend  weighting  the  data 
points  inversely  by  the  variance  of  each  measurement.  This  presumes  that  multiple 
measurements  of  the  beam  radius  are  taken  at  each  position  which  further  slows  down  the 
measurement  process.  Multi-dimensional  curve  fitting  is  not  an  easy  numerical  task  and 
available  numerical  methods  are  not  universally  robust  and  cannot  work  in  an  entirely  automated 
fashion  without  review  by  a  knowledgeable  human.  It  is  fairly  common  for  laser  technicians  to 
adjust  their  laser  or  optic  train  until  their  beam  profiler  gives  a  decent  result.  In  other  words, 
they  adjust  the  laser  to  overcome  the  problems  in  curve  fitting  and  detection  in  their  beam 
analyzer  —  opposite  from  the  ideal  case  of  adjusting  the  beam  analyzer  to  measure  the  laser  as  it 
is.  Decisions  must  be  made  on  how  to  weight  data  points,  in  what  order  to  fit  for  the  parameters, 
what  initial  guess  to  use  and  which  numerical  method  to  employ.  A  Levenberg-Marquardt^^ 
method  works  well  with  the  following  procedure:  1)  Initial  guesses  of  the  smallest  measured 
beam  radius  and  its  location  for  wo  and  zo  and  M  calculated  using  the  divergence  angle  between 
the  smallest  and  largest  measured  spot  sizes.  2)  an  unweighted  fit  on  all  three  parameters, 
keeping  the  result  for  zo  and  using  the  results  for  wo  and  as  initial  guesses  for  the  next  step. 

3)  a  fit  on  Wo  and  M  heavily  weighted  towards  points  near  focus,  keeping  the  result  for  wo  and 
using  the  result  for  as  an  initial  guess  for  the  next  step.  4)  an  unweighted  fit  on  M^.  While 
this  method  may  not  be  the  one  “best”  method,  it  is  vastly  superior  to  any  algorithm  hidden  from 
the  experimenter  by  the  word  “proprietary”.  A  known  method  can  be  analyzed  and  improved.  A 
black-box  method  will  forever  produce  results  whose  validity  is  known  only  to  the  anonymous 
developer. 

Lastly,  the  method  used  should  provide  a  means  to  assign  error  bars  based  on  signal  to  noise 
ratio  of  the  CCD  camera,  variance  in  noise  equivalent  aperture,  variance  in  beam  radius 
measurements  and  discontinuities  in  the  filter  sets.  M  is  commonly  quoted  to  2  or  3  significant 
figures  since  that  is  what  appears  on  the  displays  of  commercial  black  box  devices.  In  fact  M  is 
typically  accurate  to  slightly  better  than  1  significant  figure. 

Summary  of  noise  contributions  to  measurement  of 

Noise  on  each  pixel  of  a  CCD  camera  contributes  to  uncertainty  in  measured  beam  radius  caused 
which,  in  turn  causes  uncertainty  in  measured  M  .  The  primary  sources  of  noise  are: 


Equation  (17)  is  equivalent  to  equation  (7)  and  equation  (18)  is  equivalent  to  equation  (6)  using  the  definitions  of 
equation  (8) 
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•  Diseretization  error:  CCD  eamera  takes  diserete  measurements  of  a  eontinuous  quantity. 
Measured  as  1/resolution  of  the  eamera.  An  8-bit  eamera,  for  example,  has  a 
diseretization  error  of  1/256. 

•  CCD  noise:  Gaussian  dark  noise  present  on  all  pixels.  Measured  as  the  varianee  of 
readings  around  an  arbitrary  zero.  The  more  pixels  one  has  aeross  the  beam,  the  less  of 
an  effeet  this  eauses.  The  varianee  deereases  with  the  number  of  measurements. 

•  Filter  error:  diserete  neutral  density  filters  eause  deviation  from  saturation  of  the  camera 
and  alter  the  noise  equivalent  aperture.  Measured  as  the  minimum  percent  change  in 
filters.  If,  for  example,  the  filter  set  is  spaced  in  tenths  of  neutral  density  (ND  =.1,.2,.3,.4 
etc.),  then  this  source  of  error  is  '-'10%.  For  continuous  means  of  attenuation,  this 
contribution  is  zero. 

•  NBA  estimation  error:  NBA  (noise  equivalent  aperture)  must  be  estimated  by  using  a 
non-second  moment  method  to  measure  beam  radius.  To  the  extent  that  the  alternate 
method  differs  from  second  moment,  the  NBA  was  originally  set  incorrectly.  Due  to  dark 
current  noise  which  creates  irradiance  dependent  beam  radius  measurements,  iteration  to 
eliminate  this  source  of  error  is  not  possible.  Measured  as  a  theoretical  variance  between 
alternate  method  and  second  moment  on  a  noisy  beam. 

•  Baser  fluctuation:  Measured  as  a  variance.  Since  one  can  only  perceive  the  beam  via  a 
CCD  camera,  this  is  measured  the  same  way  as  CCD  noise,  except  with  the  laser  on. 

Total  variance  in  each  beam  radius  measurement  from  the  above  sources  is: 
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NBA=noise  equivalent  aperture 

contrast=maximum  contrast  available  on  a  particular  camera.  See  Bigure  9. 

Afilter=%  difference  between  filters. 

dark=dark  current  noise  averaged  over  the  entire  beam. 

dark-pixel=dark  current  noise  per  pixel. 

w=beam  radius. 

The  derivation  of  equation  (19)  is  presented  in  Appendix  A. 


2 

The  uncertainty  in  measuring  beam  radii  causes  an  uncertainty  in  the  M  result.  The  variance  in 

2 

M  is  equal  to  the  average  variance  in  the  data: 
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These  sources  of  error  are  far  more  significant  than  commonly  reported,  largely  due  to  the 
widespread  prevalence  of  commercial,  black  box  devices  that  include  no  estimation  of  error.  As 
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an  example,  the  author  reported  the  of  a  tunable  mid-infrared  laser  souree'^  as  “less  than  2”. 
This  was  based  on  the  above  error  analysis  and  a  statistieal  sample  of  measurements.  Based 
on  the  authors  observation  and  analysis,  measurements  more  reliable  than  1  deeimal  plaee  are 
exeeedingly  rare.  More  eommon  is  ±  0.3  for  measurements  between  1.5  and  3.  Those  who 
quote  M  =1.651,  for  example,  were  either  exeeedingly  eareful  on  an  extraordinarily  stable  laser 
or  did  no  error  analysis  at  all.  A  more  honest  number  might  be  “M^  less  than  1.8”. 

Summary  and  recommendations  for  the  use  of 

M  is  a  diffieult  measurement  to  take  properly  and  is  subjeet  to  several  theoretieal/experimental 
diseonneets.  The  aeeuraey  of  a  measurement  is  vulnerable  to  signal  to  noise  ratio  in  the  CCD 
eamera,  fluetuations  in  eontrast  ratio,  the  eontinuity  of  the  filter  set  used  and  implieations  of 
design  deeisions  regarding  eamera  resolution  and  foeusing  geometry.  M  measurements  from 
automated  “blaek  box”  instruments  are  unreliable  under  any  non-ideal  eonditions.  is 
inappropriate  for  single  shot,  multi-mode  or  annular  beams  and  is  most  appropriate  for 
laboratory,  single  mode  lasers  with  low  power  fluetuation  and  long  term  stability..  As  a  rough 
rule  of  thumb  gained  by  experienee,  M  values  above  2  are  strongly  suspeet  and  M  values  above 
5  have  lost  most  of  their  meaning.  is  best  used  below  ~1 .5  where  the  beam  has  only  one  or 
two  modes  in  more  or  less  predietable  ratios.  We  reeommend  that  a  measure  of  beam 
performanee  direetly  related  to  the  mission  objeetives  be  seleeted  rather  than  automatieally 
ehoosing  M  .  If  M  be  ehosen,  then  we  reeommend  that  identieal  measurement  apparatus  and 
proeedures  be  used  to  ensure  true  eomparability.  We  also  reeommend  that  “blaek  box”  deviees 
not  be  used  without  appropriate  error  estimation. 

-  Underlying  assumptions 

The  previous  seetion  outlined  the  experimental  diffieulties  in  measuring  M  .  This  seetion  will 
examine  the  underlying  assumptions  is  based  upon  to  enable  a  researeher  to  make  a  wiser 
deeision  regarding  the  use  of  M  for  a  partieular  appheation.  Some  measures  of  beam  quality, 
sueh  as  brightness  or  power-in-the-bueket  (depending  upon  the  way  the  ‘bueket’  size  is  ehosen), 
are  absolute  measurements,  others,  sueh  as  Strehl  ratio  and  M  are  eomparisons.  In  the  ease  of 
M^,  there  is  an  implied  eomparison  of  a  beam  with  an  embedded  beam  eomposed  entirely  of  zero 
order  Hermite-  Gaussian  profile.  Hermite-Gaussian  funetions  are  important  beeause  they  are 
self-eonsistent  solutions  of  the  paraxial  Helmholtz  equation  and  represent  shapes  whieh  will 
propagate  indefinitely.  The  mathematieal  form  of  the  magnitude  of  normalized  Hermite- 
Gaussian  funetions  is'^: 
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Where  Hn(x)  are  standard  Hermite  polynomials  and  w  (z)=Wo  (1+z  /Zr  )  with  Zr=Rayleigh 
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range  =  ttwo  IX.  The  normalized  Mn(x)  obeys  the  othogonality  relation; 
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The  orthonormal  relation  allows  the  representation  of  any  arbitrary  field  profile  E{x)  in  terms  of 
a  series  expansion  in  terms  of  these  funetions  following  the  standard  praetiees  of  linear  algebra; 

E(x)  =  Tc„u„{x,  w) 


Where  the  field  eoeffiieents  en  ean  be  determined  by  the  veetor  projeetion  of  Wn  on  E{x): 


(23) 


CO 

^  E(x)u^(x,w)dx 

—CO 


(24) 

in  one  dimension.  In  two  dimensions,  the  enm  will  be  the  produet  of  an  integral  in  x  and  another 
in  y.  Beeause  the  modes  are  orthonormal,  they  ean  be  used  as  a  basis  set  for  the  expansions  of 
equation  (23).  The  set  of  modes  Un  only  for  eomplete,  orthogonal  basis  funetions  with  respeet  to 

modes  based  upon  the  same  mode  radius,  w.  There  is  no  unique  set  of  mode  eoeffieients  that 

20 

will  deseribe  an  arbitrary  beam  shape  without  prior  determination  of  beam  radius. 

The  ‘best’  modal  decomposition 

Equation  (24)  shows  that  in  order  to  determine  the  mode  eontent  of  a  given  arbitrary  field,  one 
must  first  determine  the  eharaeteristie  beam  radius,  w.  Restrieting  our  diseussion  to  the  near 
field  at  the  plane  z=0,  that  means  we  must  first  determine  wo,  eommonly  ealled  the  beam  waist. 
The  ability  to  represent  an  arbitrary  shape  mathematieally  with  Hermite-Gaussian  modes  is 
unaffeeted  by  our  ehoiee.  There  are  many  reasonable  bases  on  whieh  to  determine  the 
eharaeteristie  beam  radius,  wo.  Eaeh  ehoiee  of  wo  will  have  a  unique  set  of  modes  assoeiated 
with  it.  All  ehoiees  of  wo  will  lead  to  a  series  representation  of  the  field.  All  will  propagate 
mathematieally  to  the  far  field.  No  physieally  measurable  properties  of  the  beam  will  be  affeeted 
by  our  ehoiee  of  Wq.  We  must  therefore  regard  the  modal  eomposition  of  a  given  beam  as  not 
entirely  a  property  of  the  beam  but  also  of  our  mathematies.  It  is  the  ease  that  wo  is  the  2"‘^ 
moment  radius  of  the  TEMoo  order  Hermite-Gaussian  mode  of  a  given  expansion.  This  will  be 
true  no  matter  what  definition  of  beam  radius  for  the  entire  beam  we  use.  The  2“‘*  moment  of  the 
lowest  mode  in  an  expansion  will  not  be  the  seeond  moment  of  the  entire  beam  exeept  in  the  ease 
of  a  pure,  TEMoo  mode,  so  a  Hermite-Gaussian  expansion  only  weakly  suggests,  but  does  not 
demand,  that  we  use  a  2"‘*  moment  definition  of  beam  radius  for  the  entire  beam.  It  remains, 
therefore,  to  make  a  good  determination  of  wo  based  on  the  physies  of  the  situation.  As  an 
illustrative  ease,  we  will  examine  the  Hermite-Gaussian  representation  of  a  beam  typieal  to  high 
energy  lasers.  Eigure  10  shows  a  square  annulus  fiat  top  field  shape,  sueh  as  one  might  obtain 
from  an  ideal  unstable  resonator.  In  this  ease  the  field  magnitude  was  arbitrarily  set  to  1,  the 
beam  has  a  width  of  5  em  and  the  inner  “hole”  has  a  width  of  2  em. 
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Figure  10  Square  Auuulus  Flat  Top 


TEMoi  +  TEMio  Hermite-Gaussian  TEM02  +  TEM20  Hermite-Gaussian  mode 

mode  (donut  mode) 

Figure  11,  Hermite-Gaussiau  modes 


Eigure  1 1  shows  several  Hermite  Gaussian  modes  which  might  be  used  as  a  basis  to  represent  the 
beam  of  Eigure  10.  The  particular  modes:  TEMoo,  TEMn,  TEM02+TEM20  and  TEMio+TEMoi 
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were  ehosen  for  this  figure  as  they  happen  to  be  the  modes  with  the  highest  field  eoeffieients  by 
several  orders  of  magnitude  for  most  ehoiees  of  oharaeteristie  beam  radius,  wo. 


Figure  12,  Field  coefficients  of  a  square  annulus  fiat  top  for  several  Hermite-Gaussian  modes  as  a  function  of 

choice  of  characteristic  beam  radius,  Wq 

Figure  12  shows  the  field  eoeffieients  for  the  modes  of  Figure  1 1  as  part  of  a  Flermite-Gaussian 
modal  deeomposition  of  the  square  annulus  flat  top  beam  of  Figure  10  as  a  funetion  of 
oharaeteristie  beam  radius.  For  referenoe,  the  vertioal  lines  in  eaoh  of  the  four  graphs  of  Figure 
12  show  the  seoond  moment  beam  radius,  whioh  is  the  ohoioe  that  M  is  based  upon.  These 
eoeffieients  were  oaloulated  for  the  z=0  plane,  but,  as  a  test,  the  square  annulus  fiat  top  was 
Fourier  propagated  to  a  distanoe  z=3000  m  and  the  same  modal  deeomposition  was  performed. 
Within  the  limits  of  numerioal  preoision,  the  shape  of  the  ourves  in  Figure  12  were  the  same  at 
z=0  and  z=3000  m.  The  point  of  this  is  to  realize  that  modal  oomposition  is  not  entirely  a 
property  of  the  beam,  but  a  produet  of  the  researeher’s  mathematies  and  ehoiee  of  oharaeteristie 
beam  radius.  Nothing  physioally  measurable  about  the  beam,  moments,  energy,  fluenoe, 
propagation  veotors,  fooal  planes  eto.  ohanges  when  one  ohooses  a  oharaeteristie  radius,  woand 
its  aooompanying  modal  deeomposition.  There  are  four  possibilities  for  determining  the  proper 
oharaeteristie  beam  radius  wo  from  the  results  shown:  first,  the  peak  of  the  ourve  of  Coo  vs.  wo; 
seoond,  the  peak  of  the  ourve  of  the  ooeffioient  oontaining  the  most  energy,  in  this  ease  C02  vs. 

Wo  (C  02+C  20  has  a  higher  peak  value  than  C  00);  third,  the  physioal  size  of  the  aperture  in  the 
near  field  and  last,  we  oould  ohoose  the  seoond-moment  or  some  other  measure  of  beam  radius  in 
the  near  field.  For  oomparison.  Figure  13  shows  the  mode  oomposition  of  a  TEMoo  (at  a 
partioular  radius)  Gaussian  beam  as  a  funetion  of  ohoioe  of  oharaeteristie  radius.  As  a  numerioal 
exeroise,  this  is  approaohing  the  absurd  sinoe  we  already  ‘know’  the  mode  oontent  of  this  beam. 
Experimentally,  however,  one  never  ‘knows’  the  mode  oontent  in  advanoe  and  must  determine  it 
from  the  data.  In  this  ease,  we  digitized  a  zero  order  Gaussian  beam  with  2"‘*  moment  radius  of 
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3.33em  from  a  hard  reetangular  aperture  20  em  wide.  Depending  upon  ehoiee  of  eharaeteristie 
beam  radius,  one  ean  have  a  signifieant  eontribution  from  the  TEM20  and  the  TEM22  mode  as 
well.  The  feature  that  truly  forees  us  to  ehoose  the  ‘eorreet’  radius  is  the  faet  that  the  Coo  curve 
reaehes  ~1.0  at  a  eertain  radius  at  whieh  all  other  modes  drop  to  zero  whieh  eorresponds  to  the 
seeond  moment  radius,  indieated  by  the  vertieal  lines.  This  is  a  luxury  we  do  not  have  with  the 
previous  example  of  Eigure  10  and  there  is  nothing  that  will  tell  us  the  ‘eorreet’  eharaeteristie 
radius,  wo,  of  a  square  annulus  flat  top. 


wO  m 

Figure  13,  Field  coefficients  of  a  TEMoo  Gaussian  beam  for  several  Hermite-Gaussian  modes  as  a  function  of 

choice  of  characteristic  beam  radius,  Wo 


and  ‘times  diffraction  limited’ 


Begin  with  equation  (23)  and  apply  the  formula  for  seeond  moment  waist,  and  use  the  identity^^ 
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Thus,  M  has  an  easy  formula  in  terms  of  Hermite-Gaussian  mode  eomposition  : 


(25) 


M^=Y,cl{2n  +  \) 


(26) 

The  last  line  of  equation  (25)  also  forms  the  general  relationship  for  the  various  “times 
diffraetion  limited”  measures.  A  measured  beam  area  in  the  far  field,  Wx  ,  is  eompared  with  an 
idealized  or  hypothesized  beam  area  whieh  is  predieted  from  a  measurement  of  beam  area  in 
the  near  field.  The  ratio  between  the  two  is  the  “times  diffraetion  limit”  of  that  beam.  This  is 
illustrated  in  the  first  line  of  equation  (27); 
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Where  W  is  the  measured  beam  radius  in  the  far  field,  rvnear  is  the  measured  beam  radius  in  the 
near  field  and  wq  is  the  ideal  waist  in  the  far  field,  z  represents  the  foeal  length  of  the  opties  and 
L  the  wavelength.  The  seeond  line  of  equation  (27)  eomes  from  referenee  [12]  and  ean  be  used 
for  illustrative  purposes  as  an  approximation  to  more  formal  means  of  mathematieally 
propagating  an  aperture  or  radius  to  a  foeal  plane.  Eaeh  ‘times  diffraetion  limited’  measure  is 
thus  tied  to  one  beam  radius  measurement  in  the  near  field  and  another  in  the  far  field.  If  the 
measure  of  beam  radius  is  ehosen  to  be  the  seeond  moment  in  both  near  and  far  field,  then  the 
‘times  diffraetion  limit’  is  equal  to  M  .  Otherwise,  it  is  not.  It  is  easy  to  imagine  that  the  ehoiee 
of  measuring  teehnique  ean  give  rise  to  a  number  of  ‘times  diffraetion  limited’  numbers  that 
eould  deseribe  a  given  beam. 


Table  3:  ‘Times  diffraction  limited’  numbers  for  square  annulus  flat-top  beam  of  Figure  10  vs.  various 
methods  of  measuring  beam  radius  in  near  and  far  field  (3km). 


far  field  i 

Coo  max® 

Largest 

mode 

maximum** 

Physical 

aperture 

2"**  moment 
radius 

R)iear  field  ^ 

3.58cm 

2.01cm 

2.5cm 

4.47cm 

Coo  max® 

2.84cm 

1 

0.32 

0.49 

1.56 

Largest 

mode 

maximum** 

5.05cm 

3.16 

1 

1.54 

5.0 

2“^*  moment 
radius 

5.37cm 

3.57 

1.12 

1.75 

5.55 

Radius  of 
central  lobe‘s 

3.03cm 

1.14 

0.36 

0.56 

1.79 

Table  3  shows  the  wide  variety  of  “times  diffraetion  limited”  numbers  we  ean  eome  up  with 
depending  upon  the  method  ehosen  to  measure  beam  radius.  We  ean  argue  that  the  square 
annulus  fiat  top  beam  of  Figure  10  is  anywhere  from  0.32  to  5.55  ‘times  diffraetion  limited’. 
How  ean  anything  be  better  than  the  diffraetion  limit?  Choose  a  suitable  basis  of  eomparison 
and  one  eertainly  ean  have  a  “better  than  diffraetion  limited”  beam!  Table  4  shows  a  similar 
analysis  for  a  zero  order  Gaussian  beam  with  a  near  field  seeond  moment  radius  of  2om. 
Through  suitable  ehoiee  of  method,  we  ean  show  that  a  zero  order  Gaussian  beam  is  from  1  to 
47.6  ‘times  diffraetion  limited’.  The  point  of  these  examples  is  not  to  argue  that  a  zero  order 


“  See  Figure  12. 

In  the  case  of  the  square  annulus  flat  top  of  Figure  10,  this  is  the  (2,0)+(0,2)  mode.  See  Figure  12. 
Measured  by  examination  of  the  far  field  pattern  after  fourier  propagation  to  the  far  field. 

©  2006  Direeted  Energy  Professional 

This  free  electronic  copy  may  be  distributed  to  individuals. 


Ross  &  Latham,  Journal  of  Directed  Energy,  2,  Summer  2006  pages,  pp.  22-58.  28/38 


beam  is  better  or  worse  than  a  square  annulus  flat  top  or  that  any  of  the  partieular  methods  of 
determining  beam  radius  ought  to  be  adopted,  but  to  emphasize  that  ‘times  diffraetion  limited’  is 
an  absolutely  meaningless  number  unless  the  methods  of  determining  beam  radius  are  defined 
along  with  it  and  that  ‘times  diffraetion  limit’  numbers  generated  by  differing  methods  of  beam 
radius  are  not  eomparable.  Differing  standards  for  measuring  beam  radius  are  not  uneommon 
and  some  are  uneonseious.  In  some  eases,  various  measures  of  beam  radius  are  identieal  for 
eertain  kinds  of  beams.  There  is  no  differenee  between  Half-Width-Half-Maximum  and  Half- 
Width- 1/e  -Maximum  for  a  flat  top  beam  in  the  aperture  plane,  but  these  measures  give  greatly 
different  results  in  the  target  plane. 

All  hard  edged  beams,  sueh  as  our  example  square  annulus  fiat  top,  have  diffraetion  ripples 
extending  to  infinity  in  the  foeal  plane.  Sinee  there  is  no  sueh  thing  as  an  infinite  deteetor,  all 
seeond  moment  measurements  on  any  hard  edged  beam  are  aetually  limited  by  the  noise 
equivalent  aperture  of  the  system,  thus  eonverting  the  measure  of  M^  to  a  elosely  related  ‘times 
diffraetion  limited’  number.  It  is  also  eommon  to  measure  far  field  beam  radius  by  width  of  the 
prineiple  lobe.  This  is  a  perfeetly  reasonable  measure  of  beam  radius,  but  eonverts  the 
measurement  from  M  to  one  of  the  other  ’times  diffraetion  limited’  numbers.  Various  beam 
radius  standards  are  listed  in  Table  3  and  Table  4  and  are  assoeiated  with  M2,  Strehl,  brightness 
and  ‘times  diffraetion  limited’.  Often,  these  measures  eonfliet  and  are  not  equivalent. 

Table  4  ‘Times  diffraction  limited’  numbers  for  a  zero  order  gaussian  beam  vs.  various  methods  of 
measuring  beam  radius  in  near  and  far  field  (3km). 


W(far 

Coo  max  or  2"‘* 

Physical 

field)! 

moment  radius 

aperture* 

fEnear  ^ 

2cm 

6cm 

Coo  max  or 

2"**  moment 

5.08cm 

1 

9 

radius 

NEA** 

1 1.68cm 

5.3 

47.6 

Summary,  Conclusion  and  Recommendations 

2 

We  have  shown  the  experimental  issues  related  to  measuring  M  along  with  the  theoretieal 
underpinnings  of  M  .  We  showed  that  the  term  “times  diffraetion  limited’’  is  meaningless 
without  referenee  to  a  partieular  method  of  determining  beam  radius  in  both  near  and  far  field. 
We  showed  that  the  theoretieal  issues  whieh  are  elear  for  nearly  zero  order  gaussian  modes 
beeome  very  eonfused  when  dealing  with  beams  typieal  of  high  power  unstable  resonators.  We 
therefore  eonelude  that  M  is  an  appropriate  measure  for  nearly  zero  order  CW  gaussian  stable 
resonator  beams  but  is  not  an  appropriate  measure  for  pulsed,  single  shot  or  unstable  resonator 
beams,  nor  ean  stable  resonator  beams  be  eompared  on  any  kind  of  equal  footing  with  unstable 
resonator  beams  in  terms  of  a  ‘times  diffraetion  limited’  number.  We  urge  that  the  term  ‘times 
the  diffraetion  limit’  be  disearded  in  any  eontext  in  whieh  laser  systems  are  eompared  against 
one  another. 


“  We  follow  the  99%  criteria  [Siegman,  Lasers,  chapter  17.1]  which  corresponds  to  a  physical  aperture  three  times  the  2"“* 
moment  radius 

Noise  Equivalent  Aperture.  See  Equation  (16)  and  Figure  9.  This  is  chosen  in  lieu  of  the  first  lobe  diameter  used  in  the 
previous  example. 
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We  also  showed  that  the  various  eommon  measures  of  beam  quality  do  not  measure  the  same 
aspeets  of  the  beam  and  that  they  should  not  beeome  eonjoined  in  our  thought  or  eommunieation. 
M^,  strehl,  ‘times  diffraetion  limit’,  brightness  and  power  in  the  bueket  do  not  measure  the  same 
thing  exeept  under  very  idealized  eonditions. 

A  mission  consistent  measure  of  optical  beam  quality 

It  is  elear  that  the  high  power  laser  eommunity  needs  a  measure  of  laser  quality  that  is  eonsistent 
with  our  mission.  An  ideal  metrie  would  satisfy  at  least  the  following  eriteria: 

•  Single  shot  measurement 

•  Taken  at  either  near  or  far  field  but  not  both 

•  Not  be  heavily  dependent  upon  noise  and  deteetor  fluetuation 

•  Not  easily  subjeet  to  obfuseation  or  argument 

•  Compare  stable  and  unstable  resonators  on  equal  footing 

•  Easily  made  uniform  for  eomparison  between  different  systems 

•  Relate  direetly  to  the  mission  requirements 

M  meets  none  of  the  above  eriteria.  Strehl  meets  several  of  the  eriteria  as  does  brightness.  The 
last  eriterion,  of  relating  to  the  mission  requirements,  depends  on  the  mission.  If  ones  mission  is 
to  have  a  very  high  peak  irradianee  without  eoneem  for  anything  else,  then  Strehl  is  appropriate. 
If  ones  mission  is  to  illuminate  a  solid  angle  (as  in  LIDAR  and  eommunieations),  then  brightness 
is  appropriate.  In  the  high  power  laser  direeted  energy  eommunity,  our  mission  is  most  often  to 
deliver  power  to  a  target.  That  is,  there  is  a  eirele  of  some  radius  and  we  wish  to  get  as  mueh 
power  in  that  eirele  as  possible  for  a  given  laser  souree.  Aside  from  the  matter  of  the  definition 
of  the  size  of  the  target  eirele,  that  is  the  definition  of  the  beam  quality  measure  known  as 
“power  in  the  bueket”  (PIB).  Even  reeognizing  the  strong  eorrelation  between  power  delivery 
applieations  and  PIB,  there  is  still  the  problem  of  the  proper  basis  of  eomparison.  PIB  is 
typieally  not  quoted  as  a  single  number,  but  shown  as  a  eurve.Eigure  14  shows  some  sample  PIB 
eurves  for  the  beams  shown  in  Eigure  15.  The  eneireled  power  (arbitrary  units)  is  shown  as  a 
funetion  of  the  far  field  radius.  In  this  ease,  all  three  eurves  were  generated  from  beams 
proeeeding  from  a  3om  diameter  aperture  without  foeusing  or  direeting  opties  to  a  distanee  of 
3km  from  the  aperture.  The  middle,  thiek  eurve  represents  the  far  field  PIB  of  an  M  =1.2  beam 
eomposed  of  0  and  E*  order  Hermite  gaussian  modes  truneated  at  the  7t/2  w  point  in  the  near 
field.  The  lower  eurve  is  for  a  donut  mode  from  the  same  aperture  and  having  the  same  energy 
eontent  as  the  low  order  Gaussian.  The  upper  eurve  is  for  a  flat  top  beam  from  the  same  size 
aperture  as  the  low  order  beam,  but  with  30%  of  its  peak  intensity  so  it  has  50%  more  energy 
eontent.  Both  stable  resonator  beams  had  woo=lcm  in  the  near  field.  A  2-dimensional 
representation  of  this  is  found  in  Eigure  11.  This  basis  of  eomparison  was  ehosen  arbitrarily.  A 
fiat  top  beam  with  the  same  peak  intensity  as  a  zero  order  Gaussian  truneated  at  the  tc/2  w  points 
has  5  times  its  energy  eontent  and  would  be  ineonvenient  to  show  on  the  same  ehart.  This  does 
emphasize  that  a  TEMoo  mode  is  very  ineffieient  at  energy  extraetion  from  most  gain  media.  If  a 
fiat  top  beam  results  in  a  loss  of  10%  diffraetion  effieieney  eompared  to  a  TEMoo  to  the  far  field 
but  delivers  5  times  more  power,  it  is  hard  to  argue  that  a  TEMoo  is  the  ‘ideal’  beam  shape. 
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Figure  14  -  PIB  curves 


comparison 


The  standard  PIB  graph  is  an  ideal  place  to  construct  a  ‘Mission  Requirements  Curve’  (MRC). 
First,  we  notice  that  lines  of  constant  average  power  run  horizontal  and  lines  of  constant  radius 
run  vertical.  If  the  physical  effect  our  mission  was  dependent  upon  were  a  function  of  average 
power  only,  we  might  stop  there.  Typically,  physical  effects  are  functions  of  intensity  or  field  on 
target.  Figure  16  shows  the  PIB  curves  of  Figure  14  with  lines  of  constant  intensity  or  field 
drawn  in.  On  a  PIB  curve,  these  lines  are  quadratic  in  shape.  For  a  given  mission,  the  minimum 
radius  is  decided  by  an  analysis  of  jitter  and  atmospheric  aberration.  The  maximum  radius  is 
decided  by  the  target  size  less  jitter.  The  physics  of  the  laser- target  interaction  determines  the 
minimum  intensity  or  field  and  a  Mission  Requirements  Curve  (MRC)  are  generated,  as  in 
Figure  17.  Note  at  this  point,  how  all  arguments  and  speculations  about  proper  basis  of 
comparison  have  been  done  away  with.  We  compare  each  beam,  not  with  a  reference  ideal,  but 
with  the  mission  requirements  based  on  the  physics  of  the  application.  Any  disputes  will  be  over 
the  mission  requirements,  not  the  beam  quality  standard. 
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Figure  16  PIB  curve  with  intensity  iso  lines 
shown 


Figure  17  PIB  curve  sample  MRC  shown 


One  unsatisfactory  element  is  left  in  the  comparison  of  a  PIB  curve  to  a  MRC  curve.  That  is,  we 
would  like  the  performance  of  the  system  to  be  reducible  to  a  single  number  so  that  one  can 
definitely  state  that  laser  system  A  is  better  than  system  B  and  by  how  much.  One  way  to  do  this 
is  to  take  the  area  of  overlap  between  a  laser’s  PIB  curve  and  the  MRC  to  generate  a  ‘Mission 
compliant  area’  (MCA).  This  is  notionally  shown  in  Figure  18.  The  two  shaded  regions  show 
the  MCA’s  for  the  flat  top  beam  and  the  Gaussian  beam.  In  the  case  chosen,  assuming  a  2kW 
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2  2 

Gaussian  beam,  the  flat  top  beam  (infinite  M  )  has  an  MCA  of  ~50  W  m  while  the  M  =1 .2  beam 
has  an  MCA  of  12.75  W  m.  In  the  ease  chosen,  we  can  conclude  that  the  fiat  top  beam  is 
approximately  four  times  a  better  beam  than  a  low  order  Gaussian  for  the  mission  envisioned 
based  on  its  overall  greater  window  of  operation  and  the  power  by  which  it  meets  or  exceeds 
mission  requirements.  The  donut  mode  beam,  in  spite  of  a  relatively  good  beam  quality  of  2.0, 
has  no  mission  compliant  area  for  this  example. 

It  will  be  observed  that  the  net  effect  of  basing  a  laser  beam  standard  on  a  previously  determined 
mission  requirements  curve  will  displace  the  center  of  argument  from  the  standard,  to  the 
mission.  Whereas  we  now  argue  about  both  the  meaning  and  propriety  of  a  laser  beam  quality 
measure  and  the  mission  requirements,  a  mission-related  standard  will  ensure  that  our 
deliberations  concern  the  mission  requirements.  The  purpose  of  a  laser  beam  quality  standard  is 
to  ensure  that  the  laser  system  will  accomplish  the  mission  if  the  standard  is  met.  Our  confusion, 
deliberations  and  debates  need  to  center  on  what  the  mission  requirements  are.  Once  that  is  over, 
the  standard  will  tell  us  if  a  given  system  will  meet  those  requirements. 
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Figure  18  Mission  compliant  areas  for  two  sample  beams 


We  realize  that  each  service  and  mission  may  have  different  standards  and  criteria  for 
constructing  a  mission  requirements  curve  and  believe  that  these  differences  are  appropriate.  It 
may  become  valuable  for  each  service  to  determine  a  small  number  of  standard  test  missions  for 
general  use  prior  to  determination  of  exact  requirements  or  for  generic  evaluation  of  new 
technologies.  We  have  only  recommended  a  general  approach  to  constructing  a  mission 
requirements  curve.  The  actual  approach,  determined  for  each  mission,  will  include  information 
from  the  laser  itself,  beam  director,  atmospheric  propagation  and  target  interaction. 


Appendix  A:  Derivation  of  error  terms  in  measurement 

Equation  (19)  represents  a  fractional  or  percent  variance  on  each  beam  radius  measurement. 
This  section  shows  the  derivation  of  its  terms  except  for  the  laser  fluctuation  and  filter  error 
which  must  be  measured  directly. 
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Discretization  and  CCD  noise: 

Take  the  seeond  moment  of  irradianee  +  some  diseretization  error  equal  to  1/resolution  of  the 
eamera.  (1/256  for  an  8-bit  eamera).  In  analyzing  the  error,  we  remember  that  the  error  oeeurs 
on  eaeh  pixel  independently  and  that  therefore  the  varianees  add  in  quadrature,  whieh  means 
they  ean  be  integrated  as  if  the  error  were  a  eonstant  under  the  integral. 
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The  limits  of  integration  are  from  the  eentroid  of  the  beam  out  to  the  noise  equivalent  aperture 
(NBA). 
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Likewise,  the  eontribution  to  beam  measurement  varianee  due  to  dark  eurrent  noise  is: 


NBA, In 


(J 


(J 


dark-beam 


X 


dark 


N,^.,NEA^ 


-rdrdO  - 


ln(eontrast)z4^,^^,o-,,,k 


ANEA^ 


The  number  of  pixels  involved  in  a  given  measurement,  Npixei=  Apixei/TrNEA^.  The  ratio  of  w  to 
NEA  is  In(oontrast).  The  error  term  in  equation  (28)  is  normalized  with  respeet  to  the  beam 
radius  and  the  radial  variable. 


NEA  estimation  error 


Uneertainty  in  the  estimation  of  the  noise  equivalent  aperture  eontributes  to  the  overall  error  in 
the  measurement  in  several  ways  in  equation  (19).  The  error  in  NEA  estimation  due  to  filter 
error,  ANEA,  is  oaleulated  as  follows,  beginning  with  equation  (16): 


NEA 

w 


.^ln(oontrast) 


NEA'  +  ANEA  =  +  .^ln(oontrast(l  +  Afilter^ )) 

^EA  =  +  Vln(l  +  Afilter") 

The  varianee  in  NEA  estimation  due  to  the  use  of  a  non-seeond  moment  method  depends,  of 
eourse,  upon  the  alternate  method  ehosen.  The  author  uses  the  full  width  at  the  1/e^  points 
(EWl/e  M).  TheEWl/e  M  is  that  point  on  a  Gaussian  where: 


©  2006  Direeted  Energy  Professional 

This  free  electronic  copy  may  be  distributed  to  individuals. 


Ross  &  Latham,  Journal  of  Directed  Energy,  2,  Summer  2006  pages,  pp.  22-58. 


33/38 


+  <^dark  +<^i 


laser 


-2 


=  e  ’ 


^method 


‘^method  ~ 


\_ 

2 

2 


1  _  ^dark  ^laser 

-1 


2  2' 
^dark  ^laser } 


Effect  of  frame  averaging  on  dark  current  noise 

It  may  be  hoped  that  averaging  over  several  shots  would  reduce  the  effect  of  dark  current  noise,  but  such  is  not  the 
case.  We  begin  by  noting  that  when  two  Gaussian  noise  distributions  are  added  together,  the  resultant  distribution  is 
the  convolution  of  the  originals.  The  convolution  of  Gaussians  of  width  g  increases  as  the  square  root  of  the  number 

of  shots  yfn  .  To  distinguish  between  averaged  distributions  and  single  shot  distributions,  the  subscript  1  will  be 
used  to  indicate  a  single  shot  quantity. 

\  C  C  ~  ^ 


Thus,  the  Odark-pkei  increases  as  the  Vn  ,  Cdark-pixeG  f^dark-pixei.i  •  This,  in  turn,  affects  the  contrast,  which  varies  as 
1  /  yfn  ,  contrast=contrastl/  -sfn  .  The  contrast  dictates  the  Noise  Equivalent  aperture  which  will  expand  as  the 

,  NEA  =  W^hl{contrast\  /  ^Jn^ .  Finally,  the  number  of  pixels  involved  in  a 
tiNEA^ 


fourth  root  of  the  number  of  shots, 
given  measurement  Npjxeis  increases  as  the  NEA^,  N 
equation  (29)and  make  the  above  substitutions: 
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The  error  in  the  2  moment  waist  caused  by  dark  current  noise  is  thus 

Aw 


w 


Apixel^ 

Icontrastl 


Whether  it  is  advisable  to  average  over  many  shots  will  depend  upon  the  noise  characteristics  of  the  laser  being 
measured.  Fortunately,  the  fourth  root  of  number  of  shots  is  a  slowly  increasing  function  and  the  benefit  of 
increasing  the  number  of  shots  is  likely  to  outweigh  the  increased  error  to  dark  current  noise  at  least  up  to  a  point. 


Sample  Calculation  of  error  in 

A  assumes  the  following^: 

•  8-bit  CCD  camera  with  the  pixel  area  of  12|umxl0|um=120|um  . 

•  contrast  =1:200. 

•  dark  current  noise  of  10%  rms,  (a^=.01). 

•  laser  noise  of  10%  rms  (a^=.01). 

•  filter  increment  of  10%. 


These  are  actual  values  from  the  Fiber  Pumped  Optical  Parametric  Oscillator  of  reference  [17]. 
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•  Assume  a  beam  radius  of  20|um  at  foeus,  for  an  NBA  of  46|um. 

•  In  the  far  field,  assume  a  1mm  beam  radius  for  an  NBA  of  2.3mm. 

Putting  these  values  into  equation  (19)  yields  a  varianee  of  (28%)  at  foeus  and  (27%)  in  the  far 
field.  The  average  varianee  will  be  approximately  (27%)^.  This  means  that  if  the  beam  analyzer 
returns  a  result  of  1.5,  we  have  a  standard  deviation  of  ±  0.4,  so  the  result  should  be  quoted  as 
1 .5±  0.4.  Alternately,  <1.9  would  be  appropriate.  In  publishing  the  results  on  the  system 
these  figures  eame  from,  we  reported  an  M  of  “less  than  2”. 

Appendix  B:  Alternate  method:  knife  edge 

Due  to  the  prevalenee  of  alternate  methods  of  M  measurement,  some  are  ineluded  in  a  brief 
appendix  of  the  ISO  standard.  This  seetion  will  rigorously  examine  the  use  of  one  of  these 
alternate  methods.  The  ISO  standard  states  that  “at  least  for  several  eases  there  exists  a 
eorrelation”  between  the  use  of  a  knife  edge  measurement  and  the  use  of  a  CCD  eamera.  The 
knife  edge  refers  to  the  praetiee  of  seanning  a  surfaee  sueh  as  a  razor  blade  aeross  the  beam,  the 
result  being  the  1-D  integral  of  the  signal  one  would  see  from  a  2-D  CCD  eamera.  The  ISO 
standard  reeommends  seanning  the  knife  until  16%  of  the  energy  reaehes  your  deteetor,  then  to 
the  84%  point  and  measuring  the  distanee  between  them.  A  faetor  is  then  multiplied  by  this 
distanee  to  ealeulate  the  “seeond  moment”  beam  radius.  Perhaps  this  method  makes  some  sense 
when  taking  measurements  by  eye  and  hand,  but  when  eonstrueting  an  automated  system  it 
raises  some  serious  issues.  Birst,  every  measurement  taken  by  an  automated  system  takes  time  as 
does  every  stage  motion.  It  may  take  a  dozen  measurements  and  a  root  finding  algorithm  to 
loeate  the  16%  point  and  the  84%  point.  Why  throw  away  all  the  extra  data?  The  seeond  issue 
is  that  “the  at  least  for  several  eases”  aetually  means  “very  nearly  single  mode  beams”. 


A  true  second  moment  knife  edge  method 

There  is  no  need  to  attempt  to  approximate  the  seeond  moment  beam  radius  using  only  two  data 
points.  Judieious  use  of  the  ealeulus  method  of  integration  by  parts  allows  the  derivation  of  a 
eonvenient  formula  that  will  use  all  the  knife  edge  data  taken; 
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In  the  ease  of  N  diserete  knife  edge  measurements,  where 
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Comparison  of  the  ISO  16-84  method  with  second  moment  knife  edge 
method 

During  the  following  discussion  we  will  refer  to  the  ISO  recommended  method  of  using  the  16% 
and  84%  points  as  the  ISO  16-84  method.  The  method  using  the  discrete  version  of  equation 
(30)  will  be  referred  to  as  the  second  moment  knife  edge  method.  We  examined  a  test  case  using 
a  modeled  beam  of  known  mode  composition  and  then  a  simulated  knife  edge  measurement 
which  will  be  processed  with  both  the  ISO  16-84  and  the  second  moment  knife  edge  methods  to 
determine  if  they  really  are  equivalent. 
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Figure  19  Comparison  of  ISO  16-84  vs.  second-moment  knife  edge 

For  the  test  case,  we  used  a  beam  composed  of  O'*’,  V\  and  4*''  order  modes.  The  can  be 
easily  calculated^^  for  comparison  by  =  Jlcliln  + 1),  where  the  Cn  are  the  energy  fractions 

of  each  mode.  Figure  19  plots  the  M  as  measured  by  the  ISO  16-84  method  and  by  the  second 
moment  knife  edge  vs.  the  calculated  for  the  test  case  mentioned.  The  second  moment 
method  used  the  numerical  curve  fitting  methods  outlined  in  the  curve  fitting  section  above.  It 
can  be  seen  that  the  second  moment  knife  edge  gives  the  known  while  the  ISO  16-84  method 
fails  to  do  so,  though  it  does  asymptotically  approach  the  proper  slope. 

Within  numerical  precision  ISO  16-84  method  consistently  under-measures  M  and  thus  fools 
the  researcher  into  thinking  he  has  a  better  beam  than  he  actually  does.  If  the  method  is  not 
specified  in  a  beam  quality  standard,  then  the  advantage  in  contract  selection  will  go  to  the 
system  that  uses  the  alternate  method  that  erroneously  reports  the  lowest  value. 

Appendix  C:  and  Strehl 

The  various  measures  of  beam  quality  all  converge  for  a  pure  zero  order  Gaussian  beam.  This 
fact  has  lead  some  to  conclude  that  the  different  measures  of  beam  quality  are  equivalent.  One 
common  version  of  this  is  to  equate  and  Strehl'^.  The  Strehl  ratio  is  the  ratio  of  the  peak 
value  of  a  beam  vs.  the  theoretical  peak  value  in  the  absence  of  any  distortions.  It  is  commonly 
used  in  imaging,  where  the  cause  of  deviation  from  zeroeth  order  Gaussian  is  due  to  atmospheric 
distortion  or  lens  aberration.  In  laser  propagation,  it  may  also  be  due  to  the  mode  content  of  the 
laser,  so  that  the  theoretical  peak  value  is  that  of  a  pure  0*’^  order  Gaussian  beam.  For  very  nearly 
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perfeet  beams,  Strehl'  is  often  approximately  equal  to  M  .  It  is  important  to  know  the 
limitations  of  this  approximation.  Both  Strehl  and  are  easy  to  ealeulate  for  a  hypothetieal 
beam  of  known  mode  eomposition.  Figure  20  shows  a  comparison  of  and  Strehl'^  for  a 
hypothetical  beam  composed  of  O*’’  order  and  1*‘  order  modes.  The  horizontal  axis  is  the  energy 
fraction  of  the  1®*  order  component  of  the  beam.  In  this  case,  and  Strehl'^  are  approximately 
equal  for  up  to  5%  1®*  order  beam,  out  to  an  of  '-1 . 1 .  Figure  21  shows  a  comparison  of 

and  Strehl'^  for  a  hypothetical  beam  composed  of  O*’’,  l’^*,  and  order  modes.  In  this  case,  the 
I  **  and  l'^  order  modes  are  given  equal  energy  content  and  the  horizontal  axis  is  the  energy 
fraction  of  higher  order  modes  in  the  beam.  In  this  case,  Strehl'^  sharply  diverges  from 
immediately  and  the  two  are  only  equivalent  for  M  =Strehr  =1.00.  In  this  particular  case,  if 
Strehl  were  measured  and  reported,  then  the  beam  quality  would  be  reported  as  better  than  it 
actually  is  and  the  laser  would  not,  in  fact,  perform  as  well  as  might  be  expected.  If  the  phrase 
“times  diffraction  limited”  were  used,  then  the  underlying  method  might  remain  obscured. 
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Figure  20  Comparison  of  and  Strehl'^  for  a  0*''  and  order  beam 
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Figure  21  Comparison  of  and  Strehl'^  for  a  O'*',  and  7*''  order  beam 
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Whether  M  is  or  is  not  approximately  equal  to  Strehl'  is  entirely  dependent  upon  the  mode 
eontent  of  the  beam.  In  praetiee,  we  ean  never  know  the  exaet  mode  eontent.  Eaeh  of  the 
different  measures  of  beam  quality  gives  a  unique  view  on  the  effeet  of  the  unknown  mode 
eontent  of  a  given  beam.  It  is  the  authors’  opinion  and  M  and  Strehf  should  never  be  assumed 
equivalent  to  one  another. 
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